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Abstract We consider the acceleration of particles due to a velocity shear in rela- 
tivistic astrophysical flows. The basic physical picture and the formation of power 
law momentum spectra is discussed for a non-relativistic velocity field using a mi- 
croscopic approach. We identify possible sites for shear acceleration in relativistic 
astrophysical jets and analyze their associated acceleration timescales. It is shown 
in particular that for a mean scattering time r scaling with the gyro-radius, the ac- 
celeration timescale for gradual shear scales in the same manner as the synchrotron 
cooling timescale so that losses may no longer be able to stop the acceleration once 
it has started to work efficiently. Finally, the possible role of shear acceleration is 
discussed with reference to the relativistic jet in the quasar 3C 273. 
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CoUimated relativistic outflows are observed in a variety of astrophysical environments, ranging 
from Active Galactic Nuclei (AGN), /i-Quasars and neutron stars in binary systems to Gamma- 
Ray Bursts (e.g., Zensus 1997; Mirabel & Rodriguez 1999; Fender et al. 2004). There is growing 
evidence today that in all these sources the relativistic outflows, or at least those parts observed 
in the radio band, are essentially launched from an accretion disk (e.g., Livio 1999; Marscher et 
al. 2002; Fender et al. 2004). In several sources, monitoring of their small-scale radio structures 
has revealed outflowing components at (apparently) high superluminal velocities suggesting 
that the bulk plasma in these jets moves at relativistic speeds and small viewing angles to- 
wards the observer. Theoretical models based on the assumption of a simple, one-dimensional 
velocity structure have thus allowed valuable insights into the emission properties of such jets. 
In general however, such models are only adequate as a first approximation, as real jets are 
likely to be characterized by a significant shear in their flow velocity field. Using theoretical and 
observational arguments relevant for AGN jets for example, we show below that at least three 
different shear flow scenarios may be distinguished. In the presence of such shear flows, particle 
acceleration is likely to occur: In their pioneer work on the particle transport in non-relativistic 
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gradual shear flows for example, Berezhko & Krymskii (1981) (cf. also Berezhko 1981) showed 
that (steady-state) power law particle momentum spectra f{p) oc p~(3+Q) j-^ay be formed if the 
scattering time follows a power law r oc p", a > 0. Later on, Earl, Jokipii & Morfill (1988) re- 
derived Parker's transport equation for non-relativistic gradual shear flows, but now augmented 
by new terms describing viscous (shear) and incrtial effects. A rclativistic generalization was 
achieved by Webb (1989) (cf. also Webb 1985 and Webb et al. 1994) based on a mixed-frame 
approach, where the scattering operator is evaluated in the comoving frame, and applied to 
the cosmic ray transport in our galaxy. Particle acceleration in non-gradual rclativistic shear 
flow has been considered by Ostrowski (1990, 1998, 2000; cf also Stawarz & Ostrowski 2002) 
based on Monte Carlo simulations, showing that very flat particle momentum spectra may be 
possible. More recently, Rieger & Mannheim (2002) have studied the acceleration of particles 
in rotating and shearing flows with application to relativistic AGN jets. Based on the insights 
gained from these works, shear acceleration in astrophysical jets appears in some respects un- 
avoidable. Whether shcar-accclcratcd particles may significantly contribute to the emission in 
a given energy band seems thus not a matter of occurrence but simply a matter of efficiency of 
the underlying shear acceleration mechanism. 

2 A SIMPLIFIED APPROACH TO SHEAR ACCELERATION 

Shear acceleration is based on the idea that energetic particles may gain energy by scatter- 
ing off systematically moving small-scale magnetic field irregularities. These scattering centres 
(SCs) are thought to be embedded in a collisionless (gradual) shear flow such that the velocity 
of a certain SC corresponds to the local flow velocity, i.e. in contrast to second order Fermi 
acceleration a random motion of SCs is neglected. Hence the acceleration process essentially 
draws on the kinetic energy of the background flow. Scattering of particles is assumed to occur 
in such a way that the particle momentum is randomized in direction and its magnitude is 
conserved in the local comoving fluid/scattering frame. Obviously, if there is no velocity shear 
(or rotation) present and the flow is not diverging, particles will neither gain energy nor mo- 
mentum merely due to scattering. However, in the presence of a velocity shear in the flow, 
the momentum of a particle travelling across the shear changes, so that a net increase may 
occur. The physical argument perhaps becomes best transparent in the microscopic picture for 
a non-relativistic continuous shear flow (cf., Jokipii & Morfill 1990) with velocity profile given 
by w = Uz{x)ez- Consider a particle with velocity vector v = {vx,Vy,Vz), relativistic mass m 
and initial momentum pi relative to local flow frame. Within a scattering time r (assumed to 
be independent of momentum) the particle will travel a distance Sx = v^t across the shear, 
so that for a gradual shear flow the flow velocity will have changed by an amount Su = Sue^, 
where 5u = (duz/dx) 6x. Contenting ourselves with a Galilean transformation for the non- 
relativistic flow speeds involved, the particle's momentum relative to the flow will have changed 
to p2 = pi + m 6u, i.e. 

P2 = pi + 2m5u pi^z + i^u)^ ■ (1) 

As the next scattering event preserves the magnitude of the particle momentum relative to 
the local flow speed, the particle magnitude will have this value in the local flow frame and 
hence a net increase in momentum may occur with time. By using spherical coordinates, defining 
< Ap/At > = 2 < (p2 — Pi) > /t as the average rate of momentum change and < Ap^ / At > = 
2 < (p2 — Pi)'^ > /t as the average rate of momentum dispersion, it can be shown that averaging 
over solid angles assuming a nearly isotropic particle distribution results in (cf., Jokipii & 
Morfifl 1990) 



Ap\ 4 fdu. 



both being related by 
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At / ~ 15 ^ V dx 



Ap\ 1 d ( 2 I A-p- 



At/ 2p^a^l^ AaT))- 



Note that a momentum-dependent scattering time obeying a power law of the form t oa 
may be accommodated by replacing 4 — > (4 + a) in Eq. (^J • We may then gain further insights 
into possible particle momentum spectra under viscous shear energy changes by studying the 
(simplified) steady state transport equation for the phase space distribution function f(p) 



p2 dp\/^ \At/ ' '^7 2p2 dp^ V \ 

assuming particles to be injected with momentum po a-nd a > 0. Solving Eq. jSJ one immediately 
obtains 

/(p)cxp-(3+") (6) 

where H{p) is the Heaviside step function. Note that for a mean scattering time scaling with 
the gyro-radius (Bohm case), i.e. t oc p, a — 1, this gives f{p) oc and thus a power law 
particle number density n{p) oc p"^ translating into a synchrotron emissivity j^, oc i>~^^'^. For 
a Kolmogorov-type {a = 1/3) or Kraichnan-type (a = 1/2) scaling on the other hand, much 
flatter spectra are obtained. 



3 POSSIBLE VELOCITY SHEAR SITES IN ASTROPHYSICAL JETS 

In general, at least three different shear scenarios may be distinguished (cf . Rieger & Duffy 2004) : 



3.1 Gradual velocity shear parallel to the jet axis 

Observationally, we have strong direct evidence for internal jet stratification (e.g. such as a fast 
moving inner spine and a slower moving boundary layer) on the pc-scale in sources such as 
3C353, M87, Mkn 501 or PKS 1055-f018 (e.g., Swain et al. 1998; Attridge et al. 1999; Perlman 
et al. 1999; Edwards et al. 2000). Moreover, there is a multitude of indirect, phenomenological 
evidence for internal jet stratification including the requirements for the unification of BL Lacs 
and FR I (e.g., Chiaberge et al. 2000) or results from hydrodynamical jet simulations (e.g., 
Aloy et al. 2000; Gomez 2002). In particular, Laing et al. (1999) have argued recently that 
the intensity and polarization systematics in kpc-scale FR I jets are suggestive of a radially 
(continuously) decreasing velocity profile Vz (r) . 

3.2 Non-gradual velocity shear parallel to the jet axis 

In general, a gradual shear analysis may be justified if the flow velocity changes continuously 
and the particle mean free path A is much smaller than the width of the shear region. However, 
when a particle becomes so energetic that its mean free path becomes larger than the width of 
the velocity transition layer, it will essentially experience a non-gradual shear flow discontinuity. 
Ostrowski (1990,1998) for example, has convincingly argued that the jet side boundary between 
the relativistic jet interior and its ambient medium may represent a natural realization of such 
a non-gradual (discontinuous) relativistic velocity shear. 
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3.3 Gradual velocity shear across the jet 

Astrophysical jets are also likely to exhibit a significant velocity shear perpendicular to their 
jet axis ('transversal shear'). In particular, several independent arguments suggest that as- 
trophysical jets may be characterized by an additional rotational velocity component: The 
strong correlation between the disk luminosity and the bulk kinetic power in the jet (e.g., 
Rawlings & Saimdcrs 1991), the successful applications of jet-disk symbiosis models (e.g., Falcke 
& Biermann 1995) and the observational evidence for a disk-origin of jets (e.g., Marscher et 
al. 2002) or periodic variability (e.g., Camenzind & Krockenberger 1992) suggest that a signifi- 
cant amount of accretion energy and hence rotational energy of the disk is channeled into the jet 
(due to angular momentum conservation). Moreover, internal jet rotation is generally implied 
in theoretical MHD models of jets as magnetized disk winds (e.g., Vlahakis & Tsinganos 1998; 
Sauty ct al. 2002). Finally, direct observational support for internal jet rotation has been re- 
cently established for stellar jets (e.g., Coffey et al. 2004). Based on such arguments, Rieger 
& Mannheim (2002) have analyzed the acceleration of particles in rotating and shearing, rela- 
tivistic AGN jets. Their analysis confirmed the results of Bcrczhko & Krymskii (1981) for the 
case of non-relativistic Keplerian rotation, but also showed that inclusion of centrifugal effects 
may lead to a flattening of the spectra. 



4 ON THE EFFICIENCY OF SHEAR ACCELERATION 

Rieger & Duffy (2004) have recently studied the acceleration potential for the different shear sce- 
narios distinguished above by estimating their associated acceleration timcscalcs. For a gradual 
shear flow parallel to the jet axis (case 3.1.) with a velocity profile decreasing linearly from rela- 
tivistic to non-relativistic speeds over a scale Ar, they found a minimum acceleration timescale 

where 7b (r) > 1 is the (position-dependent) bulk Lorentz factor of the flow and A the particle 
mean free path. By comparing the acceleration with the synchrotron cooling timescale for pa- 
rameters appropriate for pc-scale radio jets, it can be shown that efficient electron acceleration 
is quite restricted and only possible if Ar is very small, while proton acceleration is much more 
favourable since Aproion ^ Aeiectron- Perhaps most interestingly however, for a particle mean free 
path scaling with the gyro-radius, i.e. A oc 7, both the acceleration and the cooling timescale 
have the same dependency on the particle Lorentz factor 7. In such a case losses are no longer 
able to stop the acceleration process once it has started to work efficiently. Particle escape, 
when the A becomes larger than the width of the jet, or cross-field diffusion however, may then 
still represent important constraints on the maximum attainable energy. 

For a non-gradual, relativistic velocity jump (case 3.2.) on the other hand, the minimum accel- 
eration timescale is of order (cf. Ostrowski 1990, 1998) 

iacc ~ 10 — provided rg > Ar , (8) 

where denotes the gyro-radius of the particle and Ar the width of the transition layer, so 
that the acceleration timescale in this case is proportional to 7. Due to the condition Vg > 
Ar, efficient acceleration of electrons is unlikely given their associated rapid radiation losses. 
However, the acceleration of protons may be well possible until their particle mean free path 
becomes larger than the width of the jet. 

Finally, for a shear flow with relativistic Vz and azimuthal Keplerian rotation profile f2(r) oc 
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[I'ln/fY^'^ for example (case 3.3.), the acceleration timescale obeys 




(9) 



assuming the flow to be radially confined to rm<r< rj where rj is the jet radius. In general, 
efficient particle acceleration requires a region with significant rotation. For such a case the 
higher energy emission will naturally be concentrated closer to the axis (i.e. towards smaller 
radii). Again however, electron acceleration is usually very restricted (i.e. only possible for 
T ^ Hn), while proton acceleration is much more favourable (i.e. possible up to r ^ rj). If A 
scales with the gyro-radius, the acceleration and the cooling timescale will again have the same 
dependency on 7 and thus losses will not longer be able to stop the acceleration process once 
it has started to work efficiently. 



Observational and theoretical evidence suggest that astrophysical jets are characterized by a 
significant velocity shear. While internal jet rotation for example, is likely to be present at least 
in the initial parts of the jet, a signific;ant longitudinal velocity shear (parallel to the jet axis) 
prevailing all along the jet might be expected for most sources. In general, the shear acceleration 
mechanism analyzed above may thus operate over many length scales. Thus, in contrast to shock 
acceleration, shear acceleration is usually not constrained to localized regions in the jets. Besides 
the possibility that shear acceleration may account for a second population of energetic particles 
in the jet (e.g. in addition to shock-accelerated ones), shear acceleration might be of particular 
relevance for our understanding of the continuous radio and optical emission observed from 
several sources (e.g., Meisenheimer et al. 1997; Scarpa et al. 1999; Jester et al. 2001). In the 
case of the quasar 3C 273 {z = 0.158) for example, the optical spectral index is found to vary 
only very smoothly along the (large-scale) jet with no signs of strong synchrotron cooling at 
any location in the jet, e.g. between knots, contrary to expectations from shock acceleration 
scenarios (Jester et al. 2001), suggesting the need of a continuous re-acceleration mechanism 
(cf. also Meisenheimer et al. 1997). The jet in 3C 273 appears to be highly relativistic even on 
kpc-scales with a typical Doppler factor £) ~ 5 (cf. Sambruna et al 2001) so that longitudinal 
shear acceleration of particles may probably work efficiently nearly all along the jet. Moreover, 
there is also evidence for helical bulk motion in the large-scale jet (e.g., Bahcall et al. 1995) 
and internal jet helicity (e.g., oscillating ridge line, double helical pattern, periodicity) on the 
VLBI mas-scale and below (e.g., Krichbaum et al. 2000, Camenzind & Krockenberger 1992), 
suggesting that particle acceleration due to internal jet rotation (cf. Rieger & Mannheim 2002) 
may contribute on pc-scale and perhaps even on larger scales. 
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